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Abstract 

The prospects for the analysis of r channels at the SUSY benchmark point 
SPSla', especially from e+e~ — f^f^ and e~^e~ — )• t^^Tj", have been studied in full 
simulation of the ILD detector concept foreseen for the International Linear Col- 
lider. All accessible SUSY channels as well as all Standard Model backgrounds were 
generated at a centre-of-mass energy Ecms=^00 GeV, including the beam energy 
spectrum and beam backgrounds expected for nominal ILC beam parameters. With 
an integrated luminosity of 500 fb~^ delivered to the experiment, the masses of the fi 
and f2 can be measured to 107.73l[j:|]^ GeV±l.l-6M^o and 1831^^ GeV±18-(5M^o, 
respectively, while the true values in SPSla' are 107.9 GeV and 194.9 GeV, respec- 
tively. This corresponds to 5Mf^/Mf-^ ~ 0.1 % and 5Mf^/Mf^ ~ 3 % with reasonable 
assumptions on 5M^o. The cross sections for fi and f2 pair production could be 
obtained with a precision of 3.1 % and 4.2 %, respectively. Combining the mass and 
cross section measurement in the f2 channel allows to determine the LSP mass with 
a relative error of 1.7 %, assuming a known f mixing angle. In fi — )• txi decays, 
the r polarisation is measured to be 91 it 9 % and 86 it 5 % in the vr and p decay 
channels of the r, respectively. The true value in the simulated SPSla' sample is 
89.6 %. 
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1 Introduction 



The present study of the SUSY benchmark point SPSla' [T] at the International Linear 
Colhder (ILC), which was undertaken as a part of the preparation of a Letter of Intent 
for the ILD detector [2|, focuses on channels with r leptons in the final state. Contrary 
to previous fast simulation studies, e.g. |3], [1], [5], it is entirely based on full detector 
simulation of SUSY processes (signal as well as background), all Standard Model (SM) 
backgrounds, and ILC machine background. The main objective of the study was to asses 
the capabilities of the ILD detector in processes particularly sensitive to beam-beam back- 
ground and detector hermeticity. It is also a good probe for the particle identification 
capabilities and momentum resolution of the detector. Furthermore, as several fast sim- 
ulation studies have been performed in the past on similar models, the results presented 
in this work will give the opportunity to validate these results under more realistic condi- 
tions, and thus gain confidence on the validity of such studies in other channels. Finally, 
the rich phenomenology of SPSla' makes it a good model to underline the advantages of 
having an electron-positron collider that is tunable both in energy and beam polarisation. 

Specifically, the process e+e^ — t- ri^fi^ — t- t+Xi''" Xi has been studied with the goal 
to determine the achievable precision on the fi mass, the fi pair production cross section 
as well as the r polarisation in the fi decay. For the latter, the decay modes r — )■ tt^z/,- 
and T — )■ p'^Ur — TT'^TT^i'r havc been considered. In addition, the expected precision on 
the f2 mass and the f2 pair production cross section have been determined in the process 
e+e" -> r2"^f2" r+Xi^'Xi- 

The SUSY benchmark point SPSla' features a quite light mass spectrum in the slepton 
sector, and heavy squarks. Bosinos up to x-j (in e"'"e~ — ^-XiXs) would be produced at 
EcMS = 500 GeV. It is a pure mSUGRA model [6] , hence R-parity and CP are conserved. 
The unification scale parameters are: M1/2 = 250 GeV, Mq = 70 GeV, Aq = —300 GeV, 
tan/3 = 10, and sign(/i) = -|-1. The point is certainly not in contradiction with any 
experimental limits [?]. In fact, it is quite close to the most likely point indicated by 
present electro-weak precision measurements, if these are interpreted within a CMSSM 
framework [8]. In the present study, the phenomenology of SPSla' was evaluated using 
SPheno p] to run the unification scale model to the electro-weak scale. 

The fi is the next-to-lightest SUSY particle, the NLSP, with Mf, = 107.9 GeV and 
M^o = 97.7 GeV, so A(M) = 10.2 GeV. Due to this rather small mass difference the 
typical signature of e^e^ — > fi^fi^ events are two acollinear r leptons with a maximal 
energy of E^^^nax = 42.5 GeV (for Ecms = 500 GeV, and = 1.777 GeV), plus a 
significant amount of missing momentum due to the escaping neutralinos. As SPSla' 
is a point with a sizable co-annihilation contribution to the dark matter relic density, 
Mf^ is an especially important quantity to determine. This is usually done by measuring 
the upper endpoint of the energy spectrum of the r leptons from the fi decay, which is 
equal to E^-^rnax- If also the minimal allowed r energy E^-^^in can be determined, both 
Mf^ and M^o can extracted simultaneously - if not, information on M^o from other SUSY 
processes is required. In our case, £'T,mm is as low as 2.6 GeV. At such low r energies, 
the background from 77 — )• rr events is overwhelming and will pose a problem for the 
study of the fi. 
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The mass of the f2 is 194.9 GeV, so that E^-^^in = 35.0 GeV and E^- max = 152.2 GeV. 
Hence, the 77 background is less severe, but instead the background from e+e~ -^WW — )• 
iuiu becomes problematic. Another consequence of the fi being the NLSP is that r 
leptons are present in a large fraction of the SUSY decays, so that SUSY itself will be a 
major background source for r channels. 

Another important observable for characterising the f system is the r polarisation. It 
depends on both the mixing angle Of of the chiral f eigenstates into the mass eigenstates, 
as well as on the Higgsino and gaugino components of the Xi? since the interaction of 
gauginos and sfermions conserves chirality, while the Yukawa interaction of the Higgsino 
flips the chirality. The physics of the f-sector and especially of the resulting r polarisation 
is discussed in detail in [3]. 

In SPSla', the Xi is expected to have quite a small Wino component, so the x? could 
be parametrised by a single Bino-Higgsino mixing angle. To evaluate the the f mixing 
matrix, one needs to measure both f masses and 6f. This can be done by measuring 
the cross section and the endpoint of the energy spectrum of the r decay products. The 
former depends on (3^ = (1 — 4Mf^/s)^/^ and 9f, while the latter depends on Mf and M^o. 
As both f sleptons are accessible, there are four independent measurements possible to 
evaluate four parameters. However, a more sensitive channel to determine Of is the mixed 
production, e+e~ — >• fif2. To do so, it is of paramount importance to run the accelerator 
at a centre-of-mass energy between the thresholds for fif2 and f2r2 production - in SPSla' 
between 303 and 390 GeV - otherwise the background from f2f2 becomes too severe. As 
the present study is performed at Ecus = 500 GeV, we have therefore not pursued the 
determination of the mixing angle. 

The outline of the paper is as follows: We start by discussing the detector model used, 
the choice of beam polarisation, and the methods used to generate the event samples, 
followed by a breakdown of the different channels. In the following section, the method 
to extract the f signal is outlined by introducing the most significant differences between 
the signal and the various backgrounds. We continue with a description of the method 
used to determine the endpoint of the r spectrum and the cross section. Because of 
their different signatures, and different main backgrounds, the analysis for fi and f2 are 
separated into individual subsections. The subsequent section discusses the determination 
of the T polarisation, based on the r decays into vr and p, respectively. We conclude with 
a discussion on possible ameliorations to be implemented in a future study. 

2 Detector and Simulation 

The ILD detector is described in detail in [2]. For the present study the "simulation 
baseline" detector was used, which is shown in Fig. [H 

Of particular importance for the present studies are the main tracking detector (the 
TPC), the main electromagnetic calorimeter (the ECal), and the low angle calorimeters 
(the LumiCal, the LHCal and the BeamCal). 
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The track finding efficiency, even in high multiphcity events with overlaid beam back- 
ground, is 99.5 % for momenta above 1 GeV. The transverse momentum resolution 
(A(1/Pt) = A{Pt)/P^) is expected to be 2.0 x IQ-^ GeV"^ asymptotically, worsening to 
9.0 X IQ-^ GeV"^ at 10 GeV, and to 9.0 x 10"^ GeV"^ at 1 GeV. In the low angle region, 
charged tracks will be efficiently detected down to = 7°, while the only region not in the 
acceptance of the calorimetric system are the holes in the BeamCal for the beampipes. 
Around the outgoing beampipe, the radius of the hole is 20 mm at 2=3550 mm, corre- 
sponding to 5.6 mrad. Since the crossing angle of the beams is 14 mrad and the hole 
for the incoming beampipe has R=16 mm, the lower edge of the acceptance increases 
to 18.5 mrad at ~ 180°. The ECal is a highly granular SiW sampling calorimeter 
with a transverse cell size of 5 mm x 5 mm and 20 layers. In test-beam measurements 
with a prototype detector a resolution of (16.6 ± 0.1)/ a/ -E(GeV) © (1.1 ± 0.1)% has been 
achieved [2]. The simulation used here shows a similar resolution. 

The ILD_00 detector model was fully simulated using the Geant4-based Mokka [HI [12] 
package. The model not only accounts for the active elements, but also for support 
structures, for cables and cooling systems, and dead regions in the sensitive layers. In 
the TPC, the properties of the P5 gas mixture was taken into account [13], resulting in a 
parametrisation of the point error depending on both drift distance and local track-pad 
angles. Two aspects of the detector, relevant for the present study, were simulated in 
less detail: The specific energy loss in the TPC was estimated by the theoretical Bethe- 
Bloch formula, taking into account the actual TPC geometry and read-out granularity, 
to yield an estimate of the separation power between different particle species [H]. This 
estimate was used to simulate the (iii^/(ia;-based particle identification on a track- by- 
track basis. Furthermore, the response of the BeamCal to high energy electrons was 
estimated by tracking the electron to the BeamCal, and then determining the probability 
of detection from a map of the expected energy density from beamstrahlung pairs, and a 
parametrisation of the probability to detect an electron of a given energy above a given 
background energy density. Both the map and the parametrisation were obtained from a 
separate full simulation study of the BeamCal alone. 

The first means to increase the visibility of the f signal above the background is to 
determine what beam parameters are the most favourable. Because SUSY itself poses a 
background problem in the f analysis, it is required to run the ILC at the polarisation that 
minimises the background. For 100 % right e"*" polarisation and 100 % left e~ polarisation 
('P6eam(e"'", e~) = (+l,— 1)), the cross sections for X2X2 xtXi production are several 
100 fb, and the branching ratios to f are above 50 %. With the opposite polarisation, 
however, these cross sections will almost vanish. The SM cross sections are also reduced 
for Pbeam=(— 1,+1), albeit not so dramatically. In the case of the fi, an added advantage 
is that the production cross section is enhanced by a factor of 3 for Vbeam={—^,+^) 
compared to the opposite case. For the f2, the reverse is true, but the gain from the 
reduction of the background was found to be the more important feature. Hence, these 
channels were studied assuming 'P6eam=(— 0.3,+0.8), the highest level of polarisation in 
the advantageous configuration attainable with the current baseline design of the ILC. 

To generate events, the energy spectrum of the ILC beams were simulated first, in- 
cluding the effects of both the beamstrahlung and the energy spread from the main 
Linac. With these inputs, SUSY as well as Standard Model events were generated by 
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Whizard [15]. All SM channels with up to 6 fermions in the final state were simulated. 
In addition, all 8 fermions final states passing the intermediate state ttff were also gen- 
erated. For channels containing r leptons, TAUOLA (TB] was used to generate the r 
decays. By default, the helicity of the r is only internally generated in Whizard, so an 
interface between Whizard and TAUOLA was developed in order to correctly transfer this 
information between the two programs. In the case of 77 events, only multi-peripheral 
diagrams were included if the value of was below 16 GeV^. However, other types of 
diagrams (VDM, diffraction, etc.) do not produce events with large missing transverse 
momentum ^4j. Except for the channel e~^e~ — )'77ee — )• rree, a further cut at generator 
level excluded all 77 events where the invariant mass of the //-pair was below 10 GeV 
from further treatment. It should be pointed out that the simulation also includes the 
7-component of the incoming beams, so that the 77 sample also includes the processes 
when one or both of the 7's are real. 

In addition to these physics channels, the creation of e^e~ pairs due to the beam- 
strahlung process were generated, using GuineaPig [T7], and simulated with Mokka. Some 
125000 such low energetic pairs are expected to be created in each bunch crossing. A large 
fraction of these will leave the detector through the outgoing beampipe, but nevertheless 
the remaining pairs will deposit a large amount of energy in the BeamCal, and will cre- 
ate a large number of hits in the inner part of the tracking system. Approximately one 
additional charged particle per event will be detected in the TPC, and a large number of 
low energy photons will produce small showers in the calorimeters, or will convert in the 
tracking system. Due to the very large number of particles, this process cannot be simu- 
lated on an event-by-event basis, but rather a pool of 1000 bunch crossings was simulated 
separately. For each physics event, one such bunch crossing was selected at random, and 
overlaid to the event at the analysis stag^. 

The total simulated sample comprised some 13 million events. The samples gener- 
ated for the 6 and 8 fermion channels corresponded to a integrated luminosity of at least 
500 fb~^. For the four fermion channels, the sample sizes for the channels not compat- 
ible with 77 (i.e. channels not having an e"'"e~-pair in the final state) corresponded to 
between 50 and 100 fb~^, while the sample sizes of the channels compatible with 77 cor- 
responded to 0.1 fb-i (Q2 < 16 GeV^ for both beam remnants), 1 fb ^ (Q^ < 16 GeV^ 
for one beam remnant, > 16 GeV^ for the other), or 50 fb~^ {Q^ > 16 GeV^ for both 
beam remnants). Finally, for the two fermion channels, the hadronic and muonic sam- 
ples corresponded to 20 fb^^, the r+r~ sample to 500 fb~^, while the Bhabha channels 
corresponded to 1 fb^^. (In the Bhabha sample, a generator level cut of cos^^ < 0.96 and 
cos 9acoi < 0.966 was imposed). The total SM background expected at Ecus = 500 GeV 
for an integrated luminosity of 500 fb~^ was 4.92 xlO^ events, completely dominated by 
the 77 and Compton scattering processes. Apart from these two sources, the background 
was 1.53 xlO^ events, dominated by 9.94 xlO'^ of e+e^ -^ii events, mainly Bhabhas. 
The SUSY sample was divided into fifi, 7^2X2 and other channels. In all categories, the 
generated sample corresponded to at least 500 fb~^. The total number of expected events 
in the three categories were 7.92 xlO^, 8.85 xlO^, and 1.20 xlO^, respectively. 

The simulated events were reconstructed with MarlinReco [iBj. Tracks in the tracking 
system were reconstructed using the Kalman filter method, adopted from DELPHI [19j . 

^This procedure was not performed for the polarisation measurement. 
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and the calorimetric signals were treated using the particle flow paradigm, implemented 
with Pandora [20], with the PFOid package used for particle identiflcation. 

3 Mass, cross section and polarisation measurements 

The key characteristics of f production and decay, which single it out from the background, 
are: 

• only two r leptons in the flnal state 

• large missing energy and momentum 

• high acoUinearity, with little correlation to the energy of the r decay products 

• central production 

• no forward-backward asymmetry 

Different backgrounds dominate for the fi and the for the fi, the 77 background 
is important, while WW — lulu is less important; the opposite is true for the f2. For the 
mass measurement, the SUSY background is not very important, because it is dominated 
by XiXf ^-iid X2X2 production with cascade decays over f sleptons. In SPSla', the 
kinematic limits of these processes are well below those of both fi and f2 pair production, 
so they have little influence on the determination of the endpoints of the spectra. For the 
polarisation measurement, the entire spectrum is needed, and hence the SUSY background 
becomes more important. Therefore the selection cuts differ for the single analyses, but 
nevertheless two sets of common initial cuts can be deflned, one set to select the signal 
topology, and another to reject 77 events. 

3.1 Topology selection 

The flrst step in the topology selection was to group particles into jets. Jet flnding was 
performed in two ways: for the polarisation analysis, where no beam induced background 
was overlaid, the standard Durham algorithm was used [21], forced to yield two jets. In 
the mass analysis, the overlaid beam induced background implies that such a method 
has too low an efficiency due to extra tracks and clusters from the background. In this 
case, reconstructed particles to be considered for the jet clustering were flrst selected 
by demanding that all particles should have an energy exceeding 500 MeV, and that any 
charged particle should have at least one hit in the TPC To exclude neutrals reconstructed 
from showers induced by the beam background, it was demanded that they should not be 
reconstructed from hits in the BeamCal, nor should the reconstructed starting-point of the 
shower be deep in the innermost part of the forward hadronic calorimeter > 3200 mm 
and r < 500 mm). Then the algorithm used in DELPHI to flnd r leptons in f events 
was used [22] . This algorithm - to be applied only after selecting low charged multiplicity 
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events - goes through all possible combinations of charged tracks, trying to combine them 
into groups with invariant masses of Af^e* < 2 GeV, under the condition that if several 
such groupings were possible, the one with the lowest sum of jet-masses should be taken. 
In a second step, neutrals were added to the charged groups, still respecting the maximal 
mass. Any neutrals left over were labelled as belonging to the ROE group ("Rest Of 
Event"). 

Events compatible with the f topology were then selected in both analyses by requir- 
ing: 

• exactly two jets 

• less than 10 charged particle candidates 

• vanishing total charge 

• charge of each jet = ±1 

• invariant jet masses Mjet < 2.5 GeV for both jets 

• a total visible energy Ey^s < 300 GeV 

• a missing mass Mmiss > 250 GeV 

• no particle with momentum above 180 GeV 

Here Mmiss, the invariant mass of the invisible system, is determined by the difference 
between the initial e'^e^ system and the visible system. After these cuts, 9.23 xlO^ 
events remained in the standard model sample (2.41 xlO^ non-77). In the three SUSY 
subsamples, 6.92 xlO^, 7.36 xlO^ and 8.41 xlO^ events remained at this stage. 

3.2 Rejection of 77 events, and other SM processes 

The characteristics of the 77 background is the presence of two highly energetic beam 
remnant electrons at low angles, while the rest of the event forms a system of low energy 
and mass. A veto against large energy deposits at small angles was therefore imposed. 
The events that remained did so because the beam remnants escaped the detector through 
the incoming or outgoing beampipe. This limits the missing transverse momentum Pt, 
and implies that the visible system consists of two fermions almost back-to-back in the 
transverse projection. In the direction of the incoming beampipe of the other beam, a 
much larger deflection of the beam remnant could pass undetected, so a cut on Pt needs 
to depend on the azimuthal angle of the missing momentum. 

After the low angle veto, a 77 —)■ rr event could only have both high missing Pj- and 
large ^acop (the angle between the two jets projected to the plane perpendicular to the 
beam axis) if the decays of the two r leptons were asymmetric: one of the r leptons must 
have decayed with the visible products going close to opposite to the r direction - giving a 
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jet at large angle to the r, but low energy - while the other must have done the opposite - 
giving a jet with small angle to its parent r, and large energy. This property was exploited 
by the variable p_L, the scalar sum of the transverse momenta of the jets w.r.t. the thrust 
axis, in the projection perpendicular to the beam [23] . Also e^e~ -^t^t~ events tend to 
have lower p± than the SUSY signal. Figure [2] shows the simultaneous distribution of p±_ 
and $acop for signal and background events. 

This lead to the following initial cuts to reduce the 77 background: 



p_L > (2.7 sin $acop + 1-8) GeV. A similar form of this cut was proposed in [24 
no significant activity in the BeamCal 

If the azimuthal angle of the missing momentum, (pp^miss, was within 30° to the 
direction of the incoming beampipe, the value of the missing transverse momentum 
should exceed 8 GeV, see Fig. [31 



At this stage, 1.27 xlO® events remained, of which 3.82 xlO'^ were non-77 events. The 
reduction of the fifi signal sample, as well as of the SUSY background sample, is sizable 
due to these cuts: 2.26 xlO^ events remained for fifi, while 4.33 xlO^ remained for the 
background. The f2T2 sample, on the other hand, was little affected, due to its higher 
missing momentum: 6.26 xlO'^ still remained. 

In the single analyses, different cuts were used to further reduce this background to 
acceptable levels. They were intrinsically different for the various cases. The polarisation 
study relied on particle identification, and this also heavily constrained the background. 
In the f2 analysis, the signal had high enough visible energy that the 77 background could 
easily be reduced. 

In the study of the fi (both mass and polarisation in the tt channel), the following 
cuts were made in addition: 

• E^is < 120 GeV, 

• I cos 9jet\ < 0.9 for both jets, 

• $acop > 85°, 

• {Ejeti + Ejet2) sin ^acop < 30 GeV. 



The cut on {Ejeti+Ejet2) sin $acop, shown in Fig.[4|, was effective because the remaining 
SUSY background came from decays of heavier sparticles, and thus tended to be less 
back-to-back than the signal. At this stage the SM background was reduced to 9.95 xlO^ 
events, of which 1189 were non-77 events. 1.90 xlO"^ f\f\ events remained, while the 
SUSY background was reduced to 2.51 xlO^ events. 



7 



3.3 f endpoint and cross section 

The f mass can be extracted from the endpoint of the spectrum, which is equal to 
Er,max, and the Xi mass, known e.g. from a separate analysis of e and fi decays. In 
principle, the maximum of the spectrum is at ET-^min, so that the f channel can be used 
to determine M^o as well, but due to the large 77 background, the turnover point is quite 
hard to observe. 

In the analysis of both the fi and f2 mass, two additional cuts were applied against 
the 77 background: The direction of the missing momentum vector was restricted to 
I cos 6miss\ < 0.8. In events with more than 2 GeV of energy in the ROE group, at least 
20 % of that energy should be observed at angles above 30° to the beam axis. With these 
cuts, there were 2.22 xlO^ SM background events left in the sub-sample selected for the 
analysis of the fi properties (the " fi sample" in the following). At this stage it is a subset 
of the sample selected for the analysis of the f2 properties (the "f2 sample"), which still 
contains 1.15 xlO^ SM events. The remaining signal in the two cases was 1.40 xlO^ for 
fifi, and 4242 for f2f2. In addition, the fi sample contained 214 f2f2 events, and 1559 
other SUSY events, while the f"2 sample contained 1.56 xlO^ fifi events, and 3033 other 
SUSY events. 

3.3.1 fi endpoint and cross section 

For the fi sample, it was finally demanded that the mass of the visible system, M^s, 
be above 20 GeV and below 87 GeV. After these cuts, 10244 fi signal events remained 
(14.85 % efficiency), while the SM background amounts to 323.5 events. The total SUSY 
background contained 1029.6 events, including a contribution of 170.0 f2 events. Fig. [5] 
shows that the endpoint was almost background free. The turnover point, which is equal 
to the minimal r energy -E,-,mm = 2.6 GeV for most decay channels, was too distorted by 
the cuts to be measurable. 

The remaining background close to the endpoint - in the range between 30 and 70 
GeV - was found to be well described by an exponential, and was fitted in the signal-free 
region above 44 GeV. The background fit was subtracted from the observed spectrum, 
and the background-subtracted spectrum was fitted with a straight line in the range from 
30 to 41.5 GeV. As the dominating background stems from SM processes, it can be 
argued that, once ILC data is available to tune the SM simulation, the background will 
be known to a precision much better than what can be determined from the present 
simulation. Therefore, the error on the expectation value of the background was assumed 
to be negligible compared to its Poisson fiuctuations. Hence, the error on the background 
fit was not propagated into the statistical errors. Under these assumptions, the result of 
the fit was Er^max = 42.111q}2 GeV. One notes that this number is not compatible with 
the true value (42.54 GeV, with mass effects in the r decay included). The difference 
is 1.1%, equal to 3.1 a. It reflects the fact that a straight line is not quite adequate 
to describe the spectrum near the endpoint. For a flnal analysis a more sophisticated 
function would be needed, either theoretical (including the properties of the different 
decay channels and their branching ratios), or by building MC generated templates. 
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Fitting instead of the endpoint position, the result was M^^ = 107.731q;q5 GeV. 
(The true value in SPSla' is Mf^ = 107.86 GeV.) However, this result assumes that M^o 
is known. At this model point, the derivative of M^^ w.r.t. M^o is 1.1, so the actual error 
on Mf^ is 30 MeV © l.l-6M^o. In |25], an analysis of jlhp'L production using the same 
fully simulated sample was presented. The author found SM^o = 1100 MeV from this 
channel alone, so using all fi and e channels, one could expect that SM^o ^ 250 MeV. 
This is still a very conservative number, since previous fast simulation studies of /iR/iR 
and crCr production in similar, though not identical scenarios found SM^o f» 200 MeV 
and SM^o ^ 80 MeV, respectively [23]. In the most optimistic case, the relative error 
A{MfJ/Mf^ is 0.09 %, a value still dominated by 6M^o. 

The "effective" cross section, i.e. the cross section for the beam energy spectrum 
and the polarisation actually delivered by the ILC, was measured as follows: As the 
main background in this analysis arose from other SUSY channels, one could not assume 
that the background was known with arbitrarily good precision from simulation, as one 
can argue for the SM background. It needed either to be estimated from data, or to 
be reduced so much that even a very large uncertainty of its expected level had little 
influence on the final number. In principle, one could use the fact that the fi signal is 
decreased in Vteam = (+0.3, —0.8), while the SUSY background is largely enhanced to get 
an equivalent, but almost signal-free, sample. However, it was quite difficult to reduce the 
fi signal to negligible levels. This would necessitate a thorough study of the xf channel, 
which was beyond the scope of the present study. Instead, another approach was used: by 
requiring that at least one jet had an energy above 25 GeV, and none had an energy above 
42 GeV, and that no jet was identified as a single lepton, the total SUSY background 
was reduced to 49 events, accompanied by 18 SM events, while 2226 signal events still 
remained. The max(£'jet) spectrum in this sample is shown in Fig. [61 

Assuming that the uncertainty on the SUSY background is 100 %, this yielded a 
relative uncertainty on the expected number of signal events /S.{N signal) IN signal = 3.1 % 

This number of events gave an "effective" cross section aeff=Nsignai/{'C.esignai)=^58A± 
4.9 fb. The value is, by construction, equal to the expected one, and the error only comes 
from the uncertainty on N signal'- The uncertainty on C is expected to be < 0.1 % [26j, 
i.e. neghgible w.r.t. 3.1 %. For the present analysis, the error on esignai is 2 %, but this 
number could be made significantly smaller with a larger simulated sample. 

In the cross section, Mf^ only enters via an overall factor : o"^^ = A x where 

y4 is a function of Of and PbeamEl- This gives M^^ = -EfeeamA/l — (crs/yl)^/^^ and the error 
on the mass is 



^One can note that the maximum cross section occurs at Ebeam — ^5/2Mfj i.e. at 170.6 GeV. Hence, 
Ebeam= 250 GeV is not optimal for extracting the mass from the cross section: the (A^fi independent) 
1/s decrease of the cross section is more important than the rise. 
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Considering this at fixed A, 



Mf, 3(1 - 

The design goal for the ILC is < 0.1% [27], so for SPSla', and the given Ebeam, the 
second term is neghgible, and A(M^J/M^, = {A{a)/a){l3^)/3{l - 13^)= 2.1 %. 

Finally, it can be noted that, assuming the background and signal selection efficiencies 
remain the same, the optimal centre-of-mass energy for determining M^^ from the cross- 
section is 250 GeV. With these assumptions, the error on the mass is reduced to a 
quarter, i.e. A(MfJ/Mf^ = 0.5 %. It is most probable that the background levels would 
be substantially lower, since there is no SUSY background at Eqms = 250 GeV. It is 
possible that the SM background would increase, due to the change of cuts needed to 
accommodate the change in signal characteristics: -E'r,mm increases to 5.7 GeV, while 
Er^max dccrcascs to 16.9 GeV. However, this would only mean trading a poorly known 
SUSY background for a well known SM one. 

3.3.2 f2 endpoint and cross section 

For the f2 sample, it was demanded that E^s > 50 GeV, and ^acop < 155°, and the energy 
of the ROE-group was below 10 GeV. These cuts left 8095 SM background events, 3500 
of which were from WW — )■ lulu, the rest being ZZ — )• lluu or six-fermion events. In the 
f2r2 sample, 3156 events remained, while 360 SUSY background and 339 fifi events also 
passed the cuts. As the main background at this level was found to be WW — )■ lulu, it 
was requested that the other jet was not a single electron or muon, since the probability, 
in WW — lulu events, that a genuine r from one W-decaj is accompanied by an electron 
or muon on the other side is 78 %, while it is only 35 % for a r from a f decay. At this 
point, most SUSY background events had a maximum kinematically allowed jet energy 
well below the endpoint of the f2 spectrum. Only /iL/^L and clCl events where one of 
the sleptons decayed in a cascade via a Xi" to a r remained. These events were rejected 
by requiring that the most energetic jet should not be an electron or a muon. Finally, a 
likelihood ratio was formed of the joint distributions of g^eii cos9jeti and gjet2 cos9jet2 for 
signal and SM background, as shown in Fig. [3 The jet charges qjet,i are determined by the 
charge sum of the tracks in the jet, without any momentum weighting. The ratio of the 
two distributions was symmetrised along the diagonal and fitted with a multi-quadratic 
function, in order to smooth out the statistical fluctuations. The likelihood ratio exploits 
the fact that the distribution of lepton momenta in WW — )■ lulu events is forward peaked 
and forward-backward asymmetric, contrary to the signal. The final SM background was 
2257 jets in 1533 events, while the signal was 3227 jets in 1972 events (22.3 % efficiency). 
In addition 418 (4) jets in 233 (3) events survived in the fifi (SUSY background) sample. 

The endpoint of the spectrum was determined by first estimating the background. 
As virtually no SUSY background was present at the endpoint, it was assumed that 
the expected background level was known to a much better precision than its Poisson 
fluctuations. The simulated background was parametrised by an exponential function, 
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and the signal was fitted by a straight hne added to the exponential. The spectrum and 
the fitted functions are shown in Fig. [HI It should also be pointed out that only the decay 
mode f2 — )■ TXi contributes in the region where the fit was done: the cascade decays via 
a Xi^ yield jets of energies of at most 25 GeV. 

The endpoint was found to be Er,max = 150.2+^:° GeV. In this case, the simple straight 
line fit is sufficient: the value found is compatible with the true value (152.2 GeV). Fitting 
the mass gave Mf^ = 183^5^ GeV, assuming the true value for M^o. At this point, Mf^ 
is very sensitive to M^o: dM^^/dM^o = 18, so an error of 80 MeV on M^o translates into 
an additional error of 1.4 GeV on M^^- The relative error ^{Mf^) / Mf^ is 3.6 %, which is 
dominated by the error on the end-point measurement. 

To determine the "effective" cross section, a sample of events was extracted consisting 
of those where any of the two jets were accepted for the endpoint measurement. The 
events which had no jet with energy above the endpoint in fi production (42 GeV) or 
with one jet well beyond the endpoint of f2 production (160 GeV) were excluded. With 
these additional conditions, the sample contained 1765 f2 events, 1374 SM ones, and 
only 33 other SUSY events. The distribution of events is shown in Fig. |9l This yielded 
^{N signal) /N signal = 4.2 %, assumiug again a 100% uncertainty on the SUSY background, 
and a negligible uncertainty on the SM background. 

The corresponding "effective" cross section was o"e//= Nsignai/ {^^signai)= 17.7±0.7 fb, 
where once again the uncertainty on C and esignai are assumed to be negligible compared 

to A{Nsignal)- 

Also for the f2, the relative error on the beam energy is small compared to that of the 
error on the cross section, so that A(M^J/M^2 = {A{a)/a){l3^)/3{l - (3'^) = 2.4 %. The 
relative error on the mass obtained from the cross section is thus as small for f2 as for 
fi although the cross section is more than 10 times smaller and the background is much 
larger. This is due to the fact that there is almost no (poorly known) SUSY background 
in the former, and that it is on the rising edge of the cross section vs. curve at y/s = 
500 GeV. 

If, finally, the values of the endpoint and Mf^ were used to calculate M^o, one finds 
an error of 1.7 GeV, not including the error on Omix- 

3.4 r polarisation 

In this analysis, the r — )■ vr^i/,- and r — )■ p^i^r — ^ tt^tt^z/^ modes have been studied [29j. 
These modes have already been the subject of fast simulation analyses, e.g. in [3] and [lOJ. 

As explained in the introduction, the compositions of the r and the Xi manifest them- 
selves in the probability that the r has either helicity, i.e. in the r polarisation, Vr- To 
observe Vr, one uses the fact that only one helicity state exists for z/,-, which induces a 
angular distribution oc 1 — aVr cos for the visible decay products in the r rest frame [28] 
{ip is the angle between the helicity axis and the momentum, and a is a factor depend- 
ing on the mass and spin of the r decay product). The lab frame energy spectrum is 
therefore sensitive to Vt- If, in addition, the decay product is a vector, the probability of 
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the product being transversely or longitudinally polarised also depends on Vr'- one case 
would correspond to a (more likely) conservation of helicity, the other to a (less likely) 
helicity flip. Whether the vector particle is transverse or longitudinal can be determined 
by the angular distribution of its decay products: in the rest frame of a two-body decay to 
scalars, the products tend to line up along the spin axis in the longitudinal case (^ = 0), 
and to be perpendicular to it in the transversal case (^ = 7r/2). In the lab frame, the 
longitudinal case therefore corresponds to a case where the energies of the two scalars 
are maximally different, while in the transversal case they tend to be quite similar. In 
particular, the variable R = Ei/{Ei + E2) is distributed as (1 — /3yector cos{^)) , which is 
quite insensitive to the lab frame energy of the vector particle, once it is substantially 
larger than its mass (i.e. /3 ~ 1). R will be distributed as {R — 1/2)^ for the longitudinal 
case, and as 1/4 — (i? — 1/2)^ in the transversal case [3] [2H]. In SPSla', the mixing is not 
particularly larg^, and fi is expected to be mainly right-handed. Hence, is expected 
to be rather close to +1, and the spectrum in the r — )■ scalar should be harder than 
for the other helicity. In the case r — )■ vector, the vector meson is mainly longitudinal, 
yielding an R distribution peaking close to and 1. 



3.4.1 The r — > n^u^ channel 

The spectrum of the pions in the decay chain f — )■ txi — tt^i^tXi is shown in Fig. [101 
with and without ISR and beam energy spread. As the effect of these two factors clearly 
are not negligible, the true spectra were determined for extreme polarisations (-1-1 or — 1), 
and parametrised correction functions were calculated for both cases. These functions 
were double polynomials of degree 2, the two pieces being applied above or below i^T-^mm, 
respectively. With these parametrisations {F{E,+1) and F{E,—1)) at hand, the true 
spectrum for any polarisation can be obtained by applying the combined correction 

F{E, Vr) = ^-^F{E, +1) + ^-^FiE, -1) 

It should be noted that the highest sensitivity to the polarisation is in the region with 
E < E 

To extract the signal, the cuts described in Sects. 13.11 and 13.21 were supplemented 
by demanding that E^s be less than 90 GeV, and that none of the jets had an energy 
exceeding 60 GeV. The events should contain at least one signal decay candidate, defined 
as a jet that only contained a single particle, and that that particle was charged. After 
these cuts, there were 8.41x10'^ SM events (839 non 77), 201 f2f2 events, and 1678 other 
SUSY events remaining. 10730 fifi events remained, i.e. 21460 fi decays. In 4047 of 
these, the subsequent r-decay was the signal-channel r ttz/^. 

■^Due to a somewhat unlucky choice of conventions, this statement corresponds to 6mix close to 7r/2; 
The left-handed field is considered as the first chiral state, while the lighter state is considered as the 
first mass state. However, in the mass matrix, the diagonal element corresponding to the right-handed 
state is expected to be the smaller, so the lighter state in the case of the off-diagonal elements being zero 
is the pure right-handed one. Hence, with the convention, the transformation chiral state — >■ mass state 
turns the labelling "up-side down", hence a small mixing corresponds to 9 mix close to 7r/2. 
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To reject r — i- fz/^-z/^ and r — i- i^z/,- from the sample of signal candidate jets, the full 
power of particle identification of the ILD was employed: The result from the PFOid 
package, which is based on the calorimetric measurements, was supplemented by the 
measurement of dE/dx in the TPC Only about 0.4 % of the non-signal decays were 
misidentified, while the efficiency to accept signal decays was 80 %. 

This requirement was also very efficient in rejecting the remaining 77 background, 
since only a small fraction of these events did contain two r leptons. The same was true 
for 4- and 6-fermion background, albeit to a lesser extent. The background from non-fifi 
SUSY channels, on the other hand, largely contained two r leptons and was reduced only 
sightly more than the signal. It was nevertheless concluded that no further cuts were 
needed, and the final selection contained 3311 signal jets, 126 other decay modes of fi, 
334 other SUSY decays, and 122 SM jets. 

The procedure to extract the polarisation in the presence of background was to first 
fit a heuristic function to the simulated background alon^. The signal selection cuts were 
then applied to the signal+background sample, and the function describing the back- 
ground was subtracted from the observed distribution. An efficiency correction function, 
determined from signal-only simulation, was applied. As the efficiency could possibly be 
dependent on the helicities of the two r leptons in the event, the efficiency correction was 
parametrised as 



The efficiencies e (£"), e^^{E), and e^+{E), correspond to the cases of the r leptons being 

both of negative, of opposite, or both of positive helicity, respectively. These functions 
were separately determined from dedicated fast simulation samples with the corresponding 
helicity configurations. 

The ratios between initial and selected spectra are shown in Fig. [11], together with the 
fitted efficiency functions. A slight dependence on Vr was indeed observed and was found 
to be primarily caused by the cut on p±. 

The resulting distribution was then fitted with the theoretical spectrum, corrected 
for ISR and beam spread, and the polarisation was obtained, see Fig. [121 Assuming an 
integrated luminosity of 500 fb~^, the value found was Vt = (91 ± 10) %, where the error 
is statistical. The expected value in SPSla' is 89.6 %. The fitted normalisation and the 
polarisation showed a quite sizable correlation, so if the normalisation was calculated using 
the value and uncertainty of the "effective" cross section obtained in section 13.3.11 the 
error on the polarisation decreased to 6 %. The uncertainty of the average background 
stemmed from the uncertainty on its SUSY component, while the average of the SM 
component would have a negligible uncertainty. No signal-free sample with composition 
and spectrum close to that of the background in the selected sample could be constructed 

"^When real data is available, the simulation of the background can be verified by reversing cuts to 
select a signal-free, but SUSY-dominated region in the parameter space. 
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in this analysis. Hence, the background cannot be determined from the data itself, and 
one must rely on MC modelling. One could assume that no other data set than the ILC 
data will exist to validate a SUSY simulation, so the uncertainty of the model would be 
determined by the uncertainty of the ILC data itself. It is also essential that the sample 
used to verify the the SUSY simulation contains as little of the signal-channel as possible. 
Such a signal-free sample was obtained by reversing the cuts on the invariant mass of the 
other jet, the acoplanarity angle, and the cut on {Ejeti + Ejet2) sin^acop- This sample, 
shown in Fig. [121 contained 829 SUSY background jets, 128 SM jets, and 26 signal jets. 

The derivative of the fitted polarisation w.r.t. variations in the estimated SUSY 
background was determined numerically, and when multiplied by the statistical error 
on the determination of the Poisson parameter, it yielded an additional error on the 
polarisation of 5 %. 

The influence of M^^ and M^o was determined numerically, by separately varying M^o 
and -E'r.max in the fits, as the measurement of these two quantities are largely indepen- 
dent. A close to linear dependence on both these variables was found, and by using the 
uncertainty on Ermax from Sect. 13.3. H and assuming <7]\j„ ~ 250 MeV, an additional 

uncertainty of 3.4 % was determined. 
Hence, the final result was 

= 93 ± 6 ± 5 (bkg) ± 3 (SUSYmasses) % 
3.4.2 The r — )■ p^Ur — )■ n^n'^Ur channel 

In the p channel, the observable sensitive to the polarisation is E,^/Ejet which - as men- 
tioned above - is expected to be insensitive to the exact value of Ej^t, and hence to beam 
spectrum and ISR effects. Therefore, no re-evaluation of the true spectrum due to these 
effects is needed. 

Also in the p channel, the cuts described in Sects. 13. H and 13. 21 were used. In addition, 
it was demanded that Eyis be less than 90 GeV, and that none of the jets had an energy 
exceeding 43 GeV. The signal decay candidates, of which there should be at least one 
in the event, were selected by demanding that the jet only contained a single charged 
particle, and that it was accompanied by at least two neutral particles. To further reduce 
the 77 background, the cut on p± was tightened to p± > (3.5sin$acop + 2) GeV. After 
these cuts, there were 2.93x10^ SM events (733 non 77), 736 f2f2 events, and 1373 other 
SUSY events remaining. 10451 fifi events remained, i.e. 20902 fi decays. In 11120 of 
these, the following r-decay was r — > p^z/,-. 

The signal decays were selected by demanding that the corresponding jet had | cos 9 jet\ < 
0.8. The calorimeter-based PFOid algorithm was not used, because of the presence of two 
or more neutral clusters close to the track gives an unacceptably low efficiency for the sig- 
nal. The measurement of dE/dx in the TPC has no such problem, and was used to reject 

— )• e^UrT^e (accompanied by bremsstrahlung photons) and K'^Tx^Vr- Finally, the 

mass of the jet should be around the mass of the p: Mj^t G [0.4, 1.1] GeV. Only about 



14 



7 % of the non-signal decays of the Ti were misidentified, while 86 % of the signal decays 
still remained. Figure [H] shows the invariant mass spectrum of the selected events. 

The final selection contained 8165 signal jets, 1991 from other decay modes of the r in 
fifi-events (oi: 1602, K*"^: 131, all other 258), 1825 from jets in other SUSY channels, 
and 195 SM jets. In addition, the background from 77 processes was estimated to be 3000 
jets, but the lack of statistics in the simulation made it difficult to asses this number with 
precision. Due to this, the background was estimated in a somewhat less sophisticated 
manner than for the vr channel. The distribution of ET^/Ejet for the 77 before cuts was 
scaled down to correspond to the number of such events that survived all cuts, and this 
rescaled distribution was added to the background from other sources. 

Similarly to the vr channel, an efficiency correction function, determined from signal- 
only simulation, was applied. Only the efficiency is assumed be dependent on the exper- 
imental situation, not the true spectrum. Therefore, and contrary to the case of the vr 
channel, the two steps (spectrum correction and efficiency determination) could be merged 
into one, directly yielding an efficiency-corrected model prediction: 



The efficiency corrected spectra / (R), /+_(i?), and /++(i?), correspond to the cases of 

the r leptons being both of negative, of opposite, or both of positive helicity, respectively. 
These spectra were determined by fast simulation, see Fig. [151 

The fast simulation was too optimistic, both in overall selection efficiency, and the 
efficiency for low and high R. Hence the fit was restricted to R between 0.1 and 0.85, 
where the shape between full and fast simulation agreed, and the efficiency was scaled 
down equally for all polarisation configurations so that it agreed with the full simulation 
value. The observed spectrum was then fitted, with and Pr as parameters, see Fig. [161 
The result for the polarisation, simultaneously fitted with the normalisation, was found 
to be "Pt- = 86 ± 5 %. Due to the large uncertainty on the 77 contribution, it was of little 
use to study the effects of the uncertainty of the much smaller SUSY background. Due to 
the near invariance of the i?-distribution, the actual values of Mf^ and M^o are expected 
to have only a small impact on the results. 

4 Summary and Conclusions 

A study of f channels in the SPSla' SUSY scenario based on a full simulation of the ILD 
detector at the ILC was presented. 

The study was performed in the context of the detector performance studies in view 
of the ILD Letter of Intent. It was therefore based on a full detector simulation of all 
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known SM processes and machine related backgrounds. All accessible channels of the 
SPSla' SUSY model were also simulated with the same procedures. The simulation was 
done assuming that the ILC was run at a ccntre-of-mass energy of 500 GcV, delivering 
an integrated luminosity of 500 fb~^ with the electron beam being 80% right polarised, 
and the positron beam being 30% left polarised. The nominal beam parameter set was 
used to simulate the beam energy spread and beamstrahlung. 

The study has only considered f-pair production, other open channels were consid- 
ered as SUSY background. This meant that the study was done without prior knowledge 
of M^Q, so that it has concentrated on observables with low sensitivity to this param- 
eter: spectrum endpoints, cross sections and polarisation. The expected effect of the 
uncertainty on M^o on the determination of Mf is nevertheless quoted in a parametric 
form. 

Throughout, it has been assumed that the knowledge of the SM background will be 
good, so that any uncertainty on the average SM background is small compared to s/^sm 
at the final stage of event selection. The same was assumed for the determination of the 
selection efficiency. The SUSY background, on the other hand, has been assumed to be 
poorly known. It has been assigned a relative error of 100 % in most of the cases studied, 
or at best to be equal to the Poisson fluctuations in signal-free control samples, typically 
of about the same size as the final signal sample under study. 

The results on the study of the fi production for the spectrum endpoint, cross section, 
and T polarisation were: 

E,,rna. =42.11^°:^^ GcV 
^ =3.1 % 

a 

Pr =91 ±6 ±5 (bkg) ± 3 (SUSY masses) % (tt channel) 
Vr =86 ± 5 % (p channel) 

The endpoint could be used to determine Mf^, and assuming M^o has been measured to 
its nominal value (97.7 GeV) with an error of SM^o, it was found to be 

Mf^ = 107.731^:^1 ± 1-1 • SM^I GeV (endpoint). 

Also the cross section could be used to determine Mf^ . However, Ecms= 500 GeV is much 
too far from the threshold for this to be competitive: In this case A(MfJ/Mf^ would be 
2.1% assuming a known mixing angle. This should be compared to A(MfJ/Mf^ ~ 1 %o 
from the end-point, with no assumption on the mixing angle. 

The r polarisation had a lower statistical error in the p channel. However, this must be 
taken with caution, because there was a substantial amount of remaining SM background 
from 77 processes. Due to lack of simulation statistics, its contribution is poorly known. 

The results on the study of the f2 production for the spectrum endpoint and cross 
section were: 

Er,max =151.0l^;? GcV 




a 
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The endpoint value yielded 

M^, = 1831^^ ± 18 ■ 5M^i GeV (endpoint). 

For the f2, Ecms= 500 GeV is much more favourable for the determination of the mass 
from the cross section: the expected uncertainty was ^{Mr^) / Mf^ = 2.4%, comparable 
to what could be obtained from the endpoint (~ 4%). Hence, the two could be combined 
to determine M^o, and the error was found to be 1.7 GeV, similar to what was found in 
a separate analysis of /il using the same simulated sample. However, this value assumes 
that the mixing angle is known. 

The f mixing angle has not been studied in this paper, because the most sensitive 
process for its determination - fif2 production - should be studied below the f2f2 threshold 
to get a good signal to background ratio. 

In comparison with previous studies, e.g. [3], [21] and [30], several new aspects have 
been taken into account here. Most prominently, the smearing of four-vectors with design 
goal resolutions has been replaced by a detailed simulation of the various sub-detectors, 
including support structures, read-out, coolling etc. Further realism has been added by 
including not only background from Standard Model processes and from beamstrahlung 
pairs, but also from other (non-signal) SUSY processes, not always taken into account 
in the previous studies. The consideration of these additional backgrounds required im- 
provements of the r reconstruction and of the signal selection cuts. After these efforts, 
the achieved precision on the Mf^ and on the r polarisation is comparable to previous 
studies. Precise quantitative comparisons would need to take into account the different 
SUSY scenarios and accelerator parameters which have been used. The expected preci- 
sion for the f2 mass and cross section have not been evaluated in either of the previous 
studies. 

More specifically, a similar "SPSla inspired" scenario with a slightly smaller mass 
difference between fi and x? has been studied in [21] with a fast simulation of the TESLA 
detector, i.e. a predecessor to ILD. With an integrated luminosity of 200 fb~^ and a 
beam polarisation of Vbeami^^, e^) = (—0.6, -1-0.8) at a centre-of-mass energy of 400 GeV, 
statistical precisions on M^^ of 140 MeV, 100 MeV and 100 MeV have been achieved 
in the single tt, p and 37r channels, respectively, excluding any contribution from the 
uncertainty of M^o. Combined, this corresponds to a precision of about 60 MeV, quite 
similar to the Igg MeVobtained here. The higher integrated luminosity assumed in the 
study presented here is compensated in [21] by a higher degree of positron polarisation and 
more optimal choice of the centre-of-mass energy, i.e. a higher cross section. Concerning 
the r polarisation, a precision of 7 % has been achieved in [3] from the p channel. There, 
a scenario with a significantly larger mass difference has been studied using the JLC fast 
detector simulation, assuming an integrated luminosity of 100 fb~^, an electron beam 
polarisation of 0.95 and a centre-of-mass energy of 500 GeV. In view of the differences in 
beam parameters and in the SUSY scenario, this is in good agreement with the results of 
this study. 

Finally, it has to be pointed out that for many of the processes studied in this work, 
running the accelerator at Ecms= 500 GeV is not optimal. An upcoming study will treat 
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the entire SPSla' scenario as a whole, including how to partition the luminosity in an 
optimal way, and how to make use of non-f channels to measure parameters - notably 
M^o - that were found to be hard to access in the f channels. 
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Figure 1: A quadrant of the "simulation baseline" configuration of the ILD detector. 
Distances are given in mm. From [2]. 
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Figure 2: The distribution of the p± variable, defined in the text, versus the acoplanarity 
angle, a) 77 events, b) Other SM events, c) SUSY signal. The selected events are above 
the solid curve. 
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Figure 3: The distribution of the missing transverse momentum, p^missi as a function of 
the azimuthal angle of the missing momentum, (pp^miss- The accumulation of high p^miss 
events around (j)p^miss= 180° is clearly seen. It corresponds to 77 events with one beam 
remnant escaping in the incoming beampipe of the other beam. Events outside the (red) 
line were accepted. 
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Figure 4: Transverse momentum of one jet w.r.t. the direction of the other jet. Left: 
SUSY background. Right: fi signal. The cut, given by the dashed line, corresponds to 
{Ejeti + Ejf.t2) sin $acop < 30. The selected events are below the line. 
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Figure 5: The jet energy spectrum of events selected in the fi endpoint analysis, with one 
entry per jet. Light grey (yellow) histogram: signal, grey (red) : SM background, dark 
grey (green): SUSY background. The fit to the background in blue. It is fitted in the 
signal- free region to the right (solid portion of the line), and extrapolated into the signal 
region (dashed). Fit to total sample: Solid (red) line. The endpoint is at the cross-over 
of the two lines. 





Figure 6: The spectrum of the highest jet energy of events selected in the fi cross section 
analysis, with one entry per event. Light grey (yellow) histogram: signal, grey (red): 
SM background, dark grey (green): SUSY background. The cross section was estimated 
by the number of events having ma.x{Ejet) between 25 and 42 GeV, as indicated by the 
vertical lines. 
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Figure 7: Upper two figures: lego plots of qjet cos 9 j^t for the two jets, a): WW events, 
b): signal. Lower two figures: c) multi-quadratic fit to the ratio between signal and 
background. The thicker line indicates the cut: events inside this contour were accepted, 
as indicated by the arrow, d): Distribution of the likelihood ratio for signal (solid black) 
and background (dashed red). The vertical line indicates the cut. 
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Figure 8: The jet energy spectrum of events selected in the f2 endpoint analysis, with one 
entry per jet. Light grey (yellow) histogram: signal, grey (red): SM background, dark 
grey (green): SUSY background. The fit to the background is indicated by the dashed 
(blue) line. It was fitted to the background only simulation. Fit to total sample: solid 
(red) line. The endpoint is at the cross-over of the two lines. 




Figure 9: The spectrum of the maximum jet energy of events selected in the f2 cross section 
analysis, with one entry per event. Light grey (yellow) histogram: signal, grey (red): SM 
background, dark grey (green): SUSY background. The cross section was estimated by 
the number of events having ma.x{Ejet) between 42 and 160 GeV, as indicated by the 
vertical lines. 
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Figure 10: The simulated energy spectrum the pions in fi — )■ r — vr^'^z/^. The filled 
histogram shows the distribution for Eqms = 500 GeV and ISR switched off in the 
generator, while the open one shows the spectrum with the ILC beam spectrum and ISR 
included. 
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Figure 11: Ratio of the spectra before and after selection (points), and fitted efficiencies 
(lines) for different combinations of r helicity. Inverted triangles (blue): both r leptons 
are right-handed. Circles (green): the r leptons have opposite helicity. Triangles (red): 
both T leptons are left-handed. The discontinuity at 10 GeV is due to the fact that the 
particle identification algorithm changes from one set of p.d.f.s to another at that energy. 
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Figure 12: The energy spectrum of the selected pions. Light grey (yellow) histogram: 
signal. Dark grey (red) histogram: background. Black line: background fit. Grey (red) 
line: overall fit. The "dent" in the overall fit at i?,r=10 GeV is due to the discontinuity 
of the efficiency parametrisation, see Fig. [TTl 




Figure 13: Energy spectrum of the pion candidates in the control sample. Light grey 
(light green) histogram: SUSY background (non-fi). Grey (red) histogram: SM. Dark 
grey (dark green) histogram: SUSY background from other decays of the r from fi. 
Medium light grey (yellow) histogram: signal. Line: fit to spectrum. 
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Figure 14: Distribution of Mjet in the selected p sample, a): Signal and all background: 
Medium light grey (yellow) histogram: Signal. Dark grey (dark green) histogram: Back- 
ground from other decays of the r from fi. Light grey (light green) histogram: Background 
from other SUSY processes. Grey (red) histogram: Background from SM processes, ex- 
cept 77. Open dashed and solid histograms: Background from 77 — )■ rr, and 77 — )■ X. 
b): Signal and signal-type SUSY background, broken down by r decay- mode.: Light grey 
(yellow) histogram: Signal. Dark grey (green) histogram: Background from r — aiVr- 
Open histogram: Background from r — )■ Kv^-. Grey (red) histogram: Background from 
other r decays. 
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Figure 15: Distribution of i? = Et^/ Ejet before (open histogram) and after event selection 
(grey (yellow) histogram), a): both r leptons are right-handed, c): the r leptons have 
opposite helicity. r): both r leptons are left-handed. 




Figure 16: The distribution of R = Ej^/Ejet in the selected sample. The line shows the 
fitted efficiency and background corrected model. 



28 



